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ABSTRACT 



Two algorithms for finding the absolute m-center are 
developed, combining the ideas of Hakimi, Gillespie, and 
Rosenthal and Smith. The first algorithm developed is es- 
sentially a hand-computational method. It is based on 
partitioning the graph into m subgraphs centered on the 
elements of the vertex m-center. The minimum distance tree 
rooted on each element of the vertex m-center is then 
formed and modified to yield the central path and thus the 
absolute center of each subgraph. This algorithm will give 
the absolute m-centers of a graph if each of these m cen- 
tral paths passes through an element of the vertex m-cen- 
ter. The second algorithm is an iterative search of all 
possible sets of m edges on which the absolute m-center may 
be located. It is less efficient than the algorithm of 
Rosenthal and Smith when m = 1, but appears to be more ef- 
ficient for m > 1. It does eliminate the problems encoun- 
tered by the Rosenthal -Smith algorithm in handling peripheral 
vertices . 
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I. INTRODUCTION 



The optimum location problem has been examined in many 
contexts by various authors dating back as early as the 
17th Century. (See Francis [2] for an extensive biblio- 
graphy on this subject.) It arises, for example, in deter- 
mining the best location for one or more communication, 
distribution service or emergency service facilities within 
a given area to be serviced. The ’-’costs" to be minimized 
may be actual road construction costs as a function of dis- 
tance, travel time as a function of distance, the actual 
distance itself, etc. 

The location problem is modeled by a graph G(V,A): 
where V is the set of vertices which correspond to the lo- 
calities to be served, and A is the set of edges which cor- 
respond to the transportation or communication links 
interconnecting these locations. Fixed costs or distances 
are assigned to each edge of the graph. (See p. 27-30 of 
Ref. [6] for basic graph theory.) 

The m-center problem has arisen in two different forms-' 
from efforts to apply graph theoretic techniques to the 
solution of optimum location problems. The first and by 
far the most easily solved form is the vertex m-center pro- 
blem, applicable in cases where the facilities are constrained 
to be located at a vertex of the representative graph. This 
problem was defined and solved for m = 1 by Hakimi [4] , and 
was extended to cases with m > 1 by Gillespie [3] . 



5 



The research reported herein is directed toward the 
solution of the second form of the problem, the absolute 
m-center. The absolute m-center is applicable where the 
facilities to be located are only constrained to lie at 
some point of the graph G. The object then is to mini- 
mize the maximum distance from any locality to be served 
to the nearest service facility. 

This report is based on the previous works of Hakimi 
[4,5] and Gillespie [3] ; and, although conceived indepen- 
dently, it parallels some of the work of Rosenthal and 
Smith [7] . 

Section II introduces the appropriate graph- theoretic 
concepts and summarizes the solution techniques developed 
by Hakimi, Gillespie, and Rosenthal and Smith. Section III 
presents the discussion and formulation of a heuristic 
algorithm for finding the absolute center of a graph, and 
formulates but does not apply its extension to the absolute 
m-center problem. Section IV is devoted to the development 
of a second, more analytical algorithm for the solution of 
the absolute center problem and its extension to the abso- 
lute m-center. A brief example of its application is given 
in an absolute 2-center problem. Section V summarizes the 
report' and suggests some areas for further consideration. 
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II. PAST WORK 



A. THE VERTEX AND ABSOLUTE CENTERS OF A GRAPH 

Consider a connected graph G(V,A) consisting of n ver- 
tices (nodes) and M < n - 1 undirected edges (arcs) . This 
graph may be a model of a communication or transportation 
system in which traffic is allowed to flow simultaneously 
in either direction along each edge, or branch of the sys- 
tem. Let l(i,j) be the distance along the single edge (i,j) 
connecting adjacent vertices i and j, and let d(i,j) be the 
minimum distance on G between any two vertices i and j . 
Similarly, d(x,y) is the length of the shortest path on G 
between any two points x and y on G . It will also be use- 
ful to extend this notion to include several intermediate 
vertices; i.e., d(i,j,k,u,v) = d(i,j) + d(j ,k) + d(k,u) 

+ d(u,v) . Some intermediate vertices may be omitted where 
the omission will not cause confusion. 



The n x n matrix D = 



d ij 



is defined as: 



i(v i ,v .) ; i , j = 1 , . : . ,n, j f i ; 

d ij -( C1) 

^(v^v^ = 0;i = j,i = l, ... ,n. 



Define the radius associated with a vertex kcV to be 



r(k) = max d(k,v) . 
veV 

The vertex center is then defined to be that vertex in G with 
minimum radius. This radius will be called the vertex- center 



7 



radius of G, and is obtained from 



min 

r c l<j<n 



max 

l<i<n 



d C v i > V j ) 



min 

veV 



r(v) 



( 2 ) 



Similarly, the absolute center of G is defined as that point 
x* on G such that 



(3) 

min max d(v. ,x) = max d(v. ,x*) = r . 

x on G l<i<n 1 l<i<n 



1 . The Hakimi Algorithm 

In a first approach toward finding the absolute 
center of G, Hakimi [4] used equation (3) as the rationale 
for solving M simpler 'min-max' problems. The result is es- 
sentially a hand method requiring the plotting of M sets of 
linear distance functions, one set for each edge in A, to 
find M 'local centers', xi, xz, ... , x^. The value of x. 

on edge (p,q) is found by plotting the function 

(3a) 

d(v i ,x) = min{ [x + dfp^)]; [l(p,q) - x + d(q,v i )]}, 

i = 1, ... ,n, and finding the max d(v- ,x) . That x. , 

i 1 j 

j = 1, ... ,M, which minimizes (3) is then chosen as the 

absolute center, and its radius is the absolute radius. 

This method is very tedious and time consuming for 
graphs containing a large number of vertices. 

2 . ’ The Rosenthal-Smith Algorithm 

Rosenthal and Smith [6] took a very different, 
analytical approach. They began by presenting the follow- 
ing theorem, the proof of which is elementary and is omitted 
here . 
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Theorem 1 : The absolute center of graph G lies on the mid- 



point of some path which connects two vertices (not neces- 
sarily adjacent) of G. 

They define the 'central path' of G to be the 
path (v ,v ,v ) which satisfies 

A V la 



d(v x ,Vy,v z ) = min 



max d(v. ,v. ,v, ) 
i,k 1 J K 



( 4 ) 



for i,j,k = 1, ... ,n; i f k; where i and k index the columns 
and j the rows of the distance matrix D. 

The object is to find the longest path through 
each vertex of G, and then to find the vertex having the 
shortest such path. These paths must be non- recursive ; 
that is, in traveling from vertex i through j to k, each 
edge in the path is traversed only once. 

The procedure is as follows: For every row in 

D (fixed value of j) find the two largest numbers (greatest 
distances) and sum them. Call this sum d* ( v^ , Vj , v^) . Now 
find the minimum of all d* (v. , v. ,v.) as j goes from one to 

1 j K 

n. It is then necessary to check that the route found from 
i through j to k is a path for which each edge is traversed 
only once (is non-recursive). If the path found has 'back- 
paths' (is recursive), the second longest route involving j 
is found and the path checking is repeated. This process is 
continued until the longest non-recursive path through ver- 
tex j has been found. If this requirement to check each 
path for backpaths could be eliminated, the Rosenthal-Smith 
algorithm would be much more efficient. 
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The authors conclude their algorithm with the 
following theorem which describes the final step of the 
algorithm . 



Theorem 2 : The absolute center x* of graph G is located 

at the midpoint of the central path of G. 

Proof: Let the midpoint of the central path of 

G be the point x 0 . To show that Xo is the absolute center 
of G one need only show that the midpoint of any other path 
in G will yield a larger radius; because by Theorem 1 the 
absolute center lies at the midpoint of some path in G. 

The following example, originally from Hakimi 
[4] and also used by Rosenthal and Smith, will illustrate 
the latter's algorithm for the absolute center of G. 

Consider the graph in Figure 1 with the distance 
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First, compute ma> 

i ,1 

each row (vertex) . This is done by adding the two largest 
elements in each row: For row 1 this gives 24 + 34 = 48 and 

corresponds to the path (v 3 ,vi,v 5 ). These sums are listed 
in the column headed T+, shown as the first column to the 
right of the D matrix. Notice in Figure 1 that the route 
( v 3 ,vi,v 5 ) is not a non-recursive path, as the paths (v 3 ,Vi) 
and (vi ,v 5 ) have the edges (3,2) and (2,1) in common. The 



10 



next longest route involving Vi is (v 4 ,vi,Vs) with length 
54 and it is a non-recursive path; list the value 54 in 
the column headed T. The values in the T column correspond 
to the longest non-recursive path for each j . It is evi- 
dent from a comparison between elements of columns T+ and 
T that the only non-recursive path found in the first step 
(T+ column) is (vi ,vi* ,vs) . Now find the minimum entry in 
the T column. There are two rows (vertices 2 and 3) with 
the minimum value, 36; these correspond to the paths (4,2,5) 
and (4,3,1) with midpoints x and y in Figure 1, respectively. 
Both have the same absolute radius; r = r(x) = r(y) = 18. 

cl 




Figure 1: Example 1 

B. THE VERTEX M- CENTERS AND ABSOLUTE M- CENTERS OF A GRAPH 
1 . The Vertex M-Center 

Gillespie [3] extended the concept of the vertex 
center of graph G(V,A) to vertex multi - centers in the fol- 
lowing manner. A set of m vertices V * c V is called a 

m 

vertex m-center of G if, for every other set of m vertices 
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( 5 ) 



max d(v,V ) > max d(v,V *) = r ; 
veV m - veV m 



where r is defined as the radius associated with the ver- 
m 



tex m-center. 

Suppose a graph has the following distance matrix. 



D = 



0 

2 

5 

3 
5 

4 



2 

0 

3 

4 

5 

6 



5 

3 

0 

3 
2 

4 



3 

4 
3 
0 

5 
7 



5 

5 

2 

5 

0 

2 



4 

6 

4 

7 

2 

0 



Since there are six vertices there are ( 2 ) = 15 sets of two 
vertices to be examined in finding the vertex 2- center; in 
general there are (^) sets of m vertices to be examined for 
the m-center .problem. 

First, determine the minimum distance from each 
pair of vertices to every other vertex. Thus, for V 2 = {1,2} 



d ( 3 , V 2 ) = min (5,3) = 3 , 

d (4 ,V 2 ) = min (3,4) = 3 , 

d(5,V 2 ) = min (5,5) = 5 , 

d (6 , V 2 ) = min (4,6) = 4 ; 

and so on for each possible V 2 . 

Second, find the radius associated with each V 2 
where, as above, the radius is defined as the maximum of the 
set of minimum distances; i.e., the minimum distance to the 
farthest vertex from each V 2 . 

For V 2 = {1,2} the radius r(l,2) is 

r(l,2) = max d(v,V 2 ) = max (3, 3, 5, 4) = 5 . 
v 
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Third, select the set of two vertices having the 
minimum radius; i.e., 

rz = r(V 2 *) = min r(V 2 ) , (6) 

V 2 cV 

and specify that V which produces this minimum as V *. 

For the above example it may be readily verified that 
r = 3, with V 2 * = {1,5}. 

In the general m-center case equation (6) becomes 



r = mm max d(v,V ) 
m it ,, m 

V c V veV 
m 



(7) 



2 . The Absolute M-Center 

The absolute m-center of G is defined as that set 
of m points X^* on G such that for every other such set of 
m points X^ on G, 



r = max d(v,X *) < max d(v,X ) , 
am .. v ’ m J - .. v ’ nr * 

veV veV 



( 8 ) 



where r m is the absolute m-radius of G. 
am 

a. The Gillespie Algorithm 

Gillespie [3] developed an algorithm for the 
best 2-center based on the vertex 2-center and Hakimi's 
algorithm for the absolute center of a graph. He partitions 
G into- two subgraphs centered on the vertex 2- center, and 
by plotting the radius of a pair of points consisting of 
one member of the vertex 2-center and a moving point on an 
edge incident to the other member, points constituting the 
most centrally located 2-center are found. 
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Gillespie did not attempt to extend his solu- 
tion technique to cases with m > 2 , but did discuss the 
existence of the vertex and absolute m-centers and devel- 
oped the following two theorems (Theorems 1 and 2, Ref. 3 ). 

Theorem 3: Any graph containing at least m 

vertices has a vertex m-center. 

Proof: This theorem is true from 

the definition of the vertex m-center. Since, for any 

graph with m vertices there will be at least one set V , 

m’ 

there must exist a set V *, the vertex m-center. 

m 

Theorem 4: Any graph having at least m ver- 

tices has an absolute m-center. 

Proof: Theorem 4 follows directly 

from Theorem 3. Since V * always exists and V * c X , then 

m m - m 

X * must also exist, 
m 

b. The Rosenthal-Smith Algorithm 

In the extension of their algorithm for the ab- 
solute center of G to the absolute m-center, Rosenthal and 
Smith [ 7 ] partition G into m subgraphs G^ and apply their 
algorithm for the absolute center to each of these subgraphs 
An iterative procedure is then used to compare the distance 
on G from each vertex v^ , j = 1, ... ,n, to the absolute cen 

ter x^ of each G^, shifting vertices to different subgraphs 
when certain criteria are met, and recomputing the absolute 
center and radius of each affected subgraph. The algorithm 
terminates when no more vertices meet the criteria for be- 
ing shifted. 
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As the basis for their method of partitioning 
G, these authors define the 'm-node divisional path P m ' , 
where m is the number of absolute centers desired, "...as 
that path which connects m nodes (m < n) such that the dis- 
tance of the minimum branch connecting any two nodes in 
this path is greater than the minimum distance of the branch 
connecting any two nodes of any other path of G which con- 
nects m nodes." It should be observed that, as defined, 
the m-vertex (node) divisional path is actually a circuit 
when m > 2 . 

Let V m be any set of m vertices in V, then P m 
may be defined as that set V^* yielding 



max 

V eV 
m 



[ 



mm 

veV 

m 



d (v , V ) 
v ’ nr 



(9) 



where d(v.,v.) = 0 is excluded from consideration. Since 
each vertex in the path can be represented by the two dis- 
tances corresponding to the two edges of the path which are 
incident with it, (9) may be expanded in the following form 
(for four vertices) ; 



^maxj min [d(v r> v s ), d ( v s > v t )> d O t > v u )> d O u >v r )] 



The procedure for finding P m from the distance 
matrix is to find the m maximum d^ e D, i < j (upper trian- 
gular portion of D) , which form a circuit. Since all ver- 
tices of a cricuit are of degree two, no more than two entries 
may be taken from any row or column of D. The same problems 
arise here as in the procedure for finding the central path. 
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except here it must be 
deed a circuit (vice a 
once in completing the 



verified that the route found is in- 
path) , with each edge traversed only 
circuit . 
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III. THE VERTEX-M-CENTER APPROACH TO THE ABSOLUTE M-CENTER 



In most of the simple examples included in previous 
works on the m-center problem, it appears that although 
they are seldom colocated, the absolute center is frequently 
located on an edge incident to the vertex center of the 
graph. Gillespie [3] commented on this, and his efforts 
were restricted to trying to find an absolute 2-center when 
it occurs on an edge incident to a vertex 2-center. The 
algorithm developed below uses this notion as the first 
step in searching for the absolute center, but goes beyond 
it and attempts to find the absolute center even when it is 
not on an edge incident to the vertex center. 

A. DEFINITIONS AND CONVENTIONS 

The following definitions and conventions are used in 
the development of the vertex-m- center- approach algorithm: 
Consider a graph G(V,A) with vertex center c and minimum 
distance tree Tc rooted on c; i.e. , the tree such that d(v,c) 
is minimized for all v e V. Similarly, Tv will denote the 
minimum distance tree rooted on any vertex v e V. 

Define 'cross-edge' to be any edge (i,j) in A but not 
in Tc such that, if (i,j) is added to Tc, a circuit (i,c,j,i) 
is produced. 

Define 'extreme vertex' to be any vertex of degree one 
in Tc. Let the extreme vertices be ordered by distance from 
c and by the characteristics of the path to c from the ex- 
treme vertex and denoted p, q, r, ... or p^, q^, ... , 
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i = 1, 2, ... , in accordance with the following conven- 

tions : 

d(p,c) > d(q,c) > d(r,c) > ... > 0 . 

There are no edges in common between the paths (v,c), 
v = p , q , . . . ; that is , 



(p,c) H (q,c) n (r , c) ... = <*>. 



The paths (vi ,c) , (v 2 ,c), (v 3 ,c), ... have at least one 

edge in common with the path (v,c), v = p, q, r, . ..; 



(p,c) fi (p i , c) fi (p 2 ,c) ... f $ , 

(q,c) n (q i ) c) n Cq 2 ,c) ... t $ , 



and they are ordered by length 

d(p,c) > d(pi,c) > d(p 2 ,c) > ... > 0 , 

d(q,c) > d (q j , c) > d(q 2 ,c) > ... > 0 , 



Figure 2 shows a minimum distance tree rooted on the 
vertex center of a graph. Vertex 4 is the vertex center 
with a radius of ten; d(4,10) = 10 = d(p,c). The extreme 
vertices are 1, 2, 3, 5, 8, 9, and 10; with corresponding 
distances to the vertex center d(l,4) = 8, d(2,4) = 6, 
d(3 , 4) = 7, d ( 5 , 4) = 7, d(8,4) = 6, d(9,4) = 6, and d(10,4) 
= 10. Since vertex 10 is the farthest from c (vertex 4), 
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it is designated p; vertex 1 is second farthest from c and 
is designated q, etc. Vertices 3 and 5 are equidistant 
from c, but the path from 3 to c has edge (4,7) in common 
with the path from vertex 10 (p) , so vertex 3 is designated 
Pi and vertex 5 becomes r by default. There is a three-way 
tie for s between vertices 2, 8, and 9. An arbitrary choice 
is made to assign s to vertex 2 and t to vertex 9, but then 
it is necessary to assign ti to vertex 8 because its path 
to c shares the edge (4,6) with the path from vertex 9 (t) 
toe. 




Figure 2. Minimum Distance Tree Rooted on the Vertex Center 

B. AN UPPER BOUND ON THE ABSOLUTE RADIUS OF G 

Clearly the absolute radius r of G satisfies r < r , 
for if x were a local center of G such that r(x) > r c , then 
r(x) could be reduced by moving x to coincide with c, and x 
could not have been an absolute center of G. 

Consider the point x 0 on the tree in Figure 3, where c 
is the vertex center and the extreme vertices are ordered in 
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accordance with the conventions outlined in the preceding 
section. Let x 0 be the midpoint of the path (p,c,q), then 
d(x 0 ,p) = d(x 0 ,q) = r(x 0 ) = ^[d(p,c) + d(q,c)]. Since 
d(q,c) 5 d(p,c), this implies that r(x 0 ) 5 d(p,c) = r c> 
with equality holding when d(q,c) = d(p,c). It can easily 
be shown that d(x 0 ,p) > d(x 0 ,v) and d(x 0 ,p) > d(x 0 ,v^), 
v = r, s, t, ... , i = 1 , 2, ... , with equality holding 

only if d(v^,c) = d(q,c), since d(x 0 ,p) = d(x 0 ,q) by defi- 
nition. That is, Xo is no farther from any other extreme 
vertex than it is from p and q, the two most extreme ver- 
tices . 




Figure 3: A General Minimum-Distance Tree 

Now consider the point x 0 ' , Figure 3, and let d(p,x<j') > 
d(xo',a) + d(a,pi) . If it is assumed that xo ' is an absolute 
center of G (with the vertex center still at c) , then 
d(p,x 0 ') = d(q ,x 0 ’) or 

d (p jX 0 ' ) = d(x 0 ',a) + d (a , c) + d(c,q) . 

Therefore 

d(c,q) + d (c , a) < d(p,a) ; 
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but this is the condition necessary for vertex a to be the 
vertex center, which contradicts the original hypothesis 
that c is the vertex center of G. Note that the above 
analyses implicitly assume that the path (p,c,q) is the 
central path of G as defined by Rosenthal and Smith. 

Hence it is concluded that the absolute center of G 
cannot occur on an edge in Tc which is not incident to the 
vertex center if that edge is an element of the central 
path. This does not preclude its occurrence on an edge of 
an alternate tree in case of ties for vertex center, nor on 
a cross-edge excluded from Tc. 

Thus, the average of the two longest independent paths 
from the vertex center in Tc may be taken as an upper bound 
of the absolute radius of G: 

r a 5 h[d(p,c) + d (q , c) ] = r(x 0 ) . (10) 

C. MODIFYING Tc TO OBTAIN THE CENTRAL PATH 

If the path (p,c,q) on Tc (as defined above) were the 
central path of G in all cases, the problem of locating the 
absolute center of G could be solved by merely applying equa- 
tion (10) . Unfortunately this does not hold true in many 
cases. If the central path passes through a vertex center 
of G (recall that the vertex center need not be unique) it 
is frequently possible to use Tc as a starting point in 
searching for the central path. 

In general, when a cross-edge is added to Tc, alternate 
routes are formed between each pair of extremes vertices due 
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to the formation of a circuit. Thus, by inserting a cross- 
edge into Tc and removing an existing edge from Tc to break 
the circuit, a new tree Tc* is formed. If the longest path 
over Tc' is shorter than the path (p,c,q) over Tc then the 
midpoint of the newly defined path supercedes xo , the mid- 
point of (p,c,q), as a candidate for the absolute center of 
G. 

Figure 4 shows a general vertex-center tree Tc with 
representative cross-edges (shown in broken lines) which 
may be considered as alternate routes for the central path. 
The point x 0 is the midpoint of the path (p,c,q) as pre- 
viously defined. If the path (p,i,h,q) over cross-edge 
(i,h) is shorter than (p,c,q), then adding (i,h) to Tc and 
removing either edge (i,c) or (h,c) to eliminate the circuit 
(i,c,h,i) may result in a reduction of the radius of G if 
the shorter of the pair of paths, (p,i,h,c,r) and (q,h,i,c,r) 
is retained. The radius of G will be reduced if this shorter 
path is also shorter than (p,c,q). For example, assume that 
(p,i,h,c,r) is shorter than (q,h,i,c,r) and that d(p,c,q) > 
d(p,i,h,c,r) > d(p,i,h,q), and d(p,i,h,c,r) > d(q,c,r). 

Then if edge (i,h) is inserted into Tc and (i,c) is removed, 
the path (p,i,h,c,r) becomes a new candidate for the central 
path of G with a radius less than r(x 0 ) . 

Thus, if the cross-edge yielding the shortest candidate 
for the central path is found and inserted into Tc and the 
circuit is broken as outlined above, the absolute center will 
be found in those cases where the central path passes through 
a vertex- center of G. Since the vertex center need not be 



22 



unique, the minimum distance tree rooted on each of the 
alternate vertex centers must be inspected in turn. Re- 
stricting the algorithm to those cases where the central 
path does pass through a vertex center allows one to ig- 
nore the path lengths to any extreme vertices of higher 
order (closer to the vertex center) than the third most 
distant extreme vertex, r. 




Figure 4: Some Sample Cross-Edges 

To increase the generality of the discussions concern- 
ing cross-edges shown in Figure 4, it is assumed that there 
may exist additional vertices on any branch of Tc (solid 
line) in Figure 4 which are not shown in the figure. There- 
fore the distance between vertices which appear to be adjacent 
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in Tc as shown in Figure 4 will be referred to in distance 
notation rather than in edge-length notation; that is, the 
distance between vertices i and c will always be given as 
d(i,c) rather than as l(i,c) etc. There can be no inter- 
vening vertices on a cross-edge, so the length of each 
cross-edge will always be given in edge-length notation 
such as 1 (i ,h) . 

D. CROSS EDGES ON THE CENTRAL PATH 

The following paragraphs will develop the maximum al- 
lowable length for a cross-edge which will yield a reduction 
in the length of the central path if that cross-edge is in- 
serted into Tc as outlined above. This is done for each of 
the representative cross-edges shown in Figure 4. The x^, 
i = 1, ... ,8, shown in Figure 4 represent the midpoint of 

the central path candidate over the edge on which the' re- 
spective x. is located, where the candidate path is found 
as discussed in Section III.C above. Thus the point x 2 rep- 
resents the midpoint of the path (p,i,h,c,r) under the con- 
ditions of the example discussed above. Later examples will 
show that this midpoint need not fall on the cross-edge un- 
der consideration. 

First, consider the edge (k,h) with associated midpoint 
Xi in Figure 4. If (k,h) is inserted into Tc and (c,k) is 
removed, the original central path candidate (p,c,q) with 
midpoint Xo is not affected; thus no reduction in the radius 
of G is possible. If (c,h) is removed instead of (c,k) , 
then d(c,k) + l(k,h) must be shorter than d(c,h) or the path 
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(p,c,k,h,q) will not be shorter than the original path 
(p,c,q). But if d(c,k) + l(k,h) < d(c,h), the original 
path in Tc from the vertex center c to vertex h would have 
been (c,k,h) vice (c,h) by definition of a minimum distance 
tree rooted on c. This same result may be extended to edge 
(ri,qi) in Figure 4. Therefore it is concluded that no 
cross-edges between the paths (c,q) and (c,r) or (c,v), 
where d(c,v) < d(c,q) and v is an extreme vertex, need be 
cons idered . 

Now consider the cross-edge (i,h) in Figure 4. When 
(i,h) is inserted into Tc and the resulting circuit (i,c,h,i) 
is removed by breaking either of the paths (i,c) or (c,h) , 
the new path between extreme vertices p and q is (p,i,h,q). 
The length of this path must not be greater than the original 
path between p and q, (p,c,q), or the radius of G will be in- 
creased; and it must be less than (p,c,q) to reduce the 
radius of G. In addition, the longer of the two paths 
(p,i,h,c,r) and (q,h,i,c,r) can be removed by breaking either 
(c,h) or (c,i) respectively. Therefore the shorter of these 
two paths will remain and must also be shorter than (p,c,q) 
to allow a reduction in the radius of G. These conditions 
are equivalent to 

l(i,h) < d(i,c,h) = d(i,c) + d(c,h) (11) 

and 



d(p,i) + 1 (i ,h) + d (h , c , r)) 



mm < 



, d(q,h) + 1 (i ,h) + d(i,c,r), 



( 12 ) 
5 
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Inequality (12) may be rewritten in the same form as (11) , 



f d (i , c) + d (h , q) - d(c,r) 

1 (i ,h) < max s . (12a) 

l d (h , c) + d(i ,p) - d(c,r) 



Therefore, if l(i,h) satisfies inequalities (11) and (12), 
cross-edge (i,h) may be inserted into Tc, the resulting 
circuit removed by breaking the longer of the two paths 
(p,i,h,c,r) and (q,h,i,c,r) between c and i or h, and a new 
candidate for the central path will be obtained. If strict 
inequality holds in both (11) and (12) , this will yield a 
new radius of G which is strictly less than r(x 0 ) . 

Next, consider the cross-edge (j ,h) in Figure 4; (j ,h) 
differs from (i,h) in that the secondary extreme vertex p x 
is connected to the path (p,c) at a point (vertex e) which 
is closer to the vertex center than is vertex j . It is 
therefore possible for the path (pi ,e ,j ,h,c,r) to be longer 
than the path (p,j,h,c,r), and inequality (12) must be modi- 
fied to hold for the longer of these two paths. When this 
is done, (12) becomes 



1 (i ,h) < max 



j d ( p , j , e , c) + d (h , q) - d ( c , r )' 
<1 -max [d(px ,e , j) ; d(p,j)] _ 



(12b) 



^ d(p , j ) + d (h , c) - d(c,r) 



Inequality (11) remains in the same form, with j sub- 



stituted for i; 



l(j,h) < d ( j , c ,h) = d ( j , c) + d(c,h) (Ha) 



The roles of p and p x may be interchanged in this case, 
and similar branching may occur in paths associated with 
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with vertex £ such as over cross-edges (i,f) or (j,f) (not 
shown in Figure 4) . These situations and role interchanges 
between q and qi would still be handled in a similar manner 
but modifying (11) and (12) as was done in (11a) and (12b) . 

The same procedures also apply for the cross-edge (j ,k) 
in Figure 4. In this case the roles of q and r are inter- 
changed in (12a), and (11) is modified also. 

To consider the cross-edges (p,q) and (p,r) in Fig- 
ure 4, the simplified graph in Figure 5 will be useful. 

For cross-edge (p,q) , in the absence of secondary extreme 
vertices as shown in Figure 5, the appropriate form of the 
path length criterion (12) becomes 

T l(P,q) + d(q,c) + d ( c , r )*! 

(a) min k< d(p,c) + d(q,c) , 

[ 1 (p , q) + d(p,c) + d(c,r)J 

and the appropriate form of (11) is 

(b) l(p,q) 5 d(p,c) + d(q,c) . 

But (b) is less restrictive than (a) , and so can be discarded. 
Then inequality (a) may be re-written in the form 

f d(P>c) 

(ai) 1 (q ,p) + d(c , r) < max < 

l d(q>c) 

But by 'definition , d(p,c) > d(q,c) ; therefore (ai) becomes 

l(p,q) < d (p , c) - d(c,r) . (12c) 

For the cross-edge (p,r), in the absence of secondary 
extreme vertices as in Figure 5, the appropriate form of path 
length criterion (12) is 
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min 



(b) 



l(p,r) + d(c,r) 



+ 



d ( c ,q)' 

► 



< 



d(c,p) 



+ d(c,q) 



v 1 (p , r) + d(c,p) + d(c,q)J 

and since d(c,p) - d(c,r) > d(c,p) - d(c,q) by definition, 
this reduces to 



1 (P > r) 5 d(c,p) - d(c,r) . (12d) 

Notice that the right sides of (12c) and (12d) are identical 
This is due to the simple role interchange of q and r in the 
left side of (a) to arrive at (b) , and illustrates the ease 
with which the situations discussed on the preceding page 
may be handled. Once again the proper form of (11) is less 
restrictive than (b) and may be ignored, but in this case 
(12d) is^ the appropriate form of (11) . When secondary ver- 
tices are involved, the appropriate form of (12) is not so 
simple as (12c) and (12d) , as is shown below. 




Figure 5: Simplification of Figure 4. 

Now return to Figure 4 and reconsider cross-edge (p,q) . 
In the presence of secondary extreme vertices, the proper 
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form of distance criterion (12) is 



(c) 



min< 



d (Pi>e,p) + l(p,q) + d(q,c,r)' 
d (q ! ,b ,q) + l(p,q) + d(p,c,r)J 



5 d (P',c,q) ; 



If d(pi,e,p) = 0 and d(qi,b,q) = 0, (c) reduces to (12c) 
as would be expected. As in the previous case, (c) is more 
restrictive than the appropriate form of (11) , and is 
equivalent to 



d(e,c) - d(pi ,e) 



l(p,q) 5 max < 



- d(c,r) . (12e) 



d(b,c) - d (q i ,b) 

From (c) one may also deduce the further dual requirements 
that 

d(p,c) > d(p ! ,c) , 



(13) 



d (q >c) > d (q i ,c) 



By interchanging the roles of p and px and of q and q x , 
inequality (c) becomes applicable to a cross-edge between 
secondary vertices such as (pi,qi) in Figure 4. In this 
case the proper form of (11) is no longer less restrictive 
than the revised version of (c) , and so must be considered 
again. ' Relation (11) now becomes 



(lib) 

l(Pi,qi) < [ d ( e ,c) - d (e ,p x ) ] + [d(b,c) - d(b, qi )] 
While the correct form of (c) is 

f [d(e , c) - d(e ,p i ) ] + [d (b , q) - d(b , qi ) ]) (12£) 

< (" d ( c , r 

[d(b,c) - d(b ,qi) ] + [d(e ,p) - d(e, P x)]) 



1 (P l >qi) f max 
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Throughout the preceding discussion, one important fac- 
tor has been ignored, primarily because it is impossible to 
quantify in a general sense. When a cross-edge has been ad- 
ded to Tc and an edge of Tc is removed to break the resul- 
tant circuit, additional extreme vertices which were not 
extreme in Tc may be generated in the modified tree. For 
instance, if cross-edge (j ,h) in Figure 4 has been found to 
satisfy (11a) and (12b) and is added to the tree in Figure 4, 
when the circuit (j ,e,g,i ,c,h, j) is broken by removing (g,i) 
(assuming path (g,i) to now be a single edge), vertices g 
and i become newly defined extreme vertices. It is then pos- 
sible for the path (i,c,h,j,p) or (g ,e, j ,h,c,r) to be longer 
than the initial controlling path, (p,c,q). This possibility 
must therefore be avoided in all instances. Example 4, 
Section III.F will illustrate a case in which it is not pos- 
sible to modify Tc because of the occurrence of newly de- 
fined extreme vertices with associated path lengths greater 
than the length of (p,c,q). 

E. THE ALGORITHM 

The discussion of the preceding sections may be sum- 
marized in step-wise manner as follows: 

1.. Find the vertex center of the graph. See Refs. 

[3] and [4] or Section II above. 

2. Form the minimal distance tree Tc rooted on the 
vertex center c using the Dijkstra Algorithm or a similar 
method (see Dreyfus [1]) . 
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3. Order the extreme vertices as outlined in Section 
III. A. 

4. Find the midpoint Xq of the longest path (p,c,q) 
over Tc. 

5. Check all cross-edges from vertices in the path 
from p to c with the applicable form of inequalities (11) 
and (12) . 

6. Choose the cross-edge which yields the shortest 
central path. 

7. Insert that cross-edge into Tc and break the re- 
sulting circuit so that the longest non- recursive path be- 
tween a pair of extreme vertices is removed. Break the 
circuit so as to minimize the length of the longest such 
path from a newly defined extreme vertex. 

8. If the longest path from a newly defined extreme 
vertex is longer than path (p,c,q), drop this cross-edge 
from consideration, return to step 6 to find the cross-edge 
yielding the next shortest central path candidate. 

9. If the longest path from a newly defined extreme 
vertex is shorter than (p,c,q) but longer than the longest 
path defined by the applicable forms of (11) and (12) , re- 
tain the result but return to 6 to find the cross-edge 
yielding the next longer central path candidate. 

10. Find the midpoint of the shortest central path 
candidate (longest remaining path between any pair of ex- 
treme vertices including newly defined extreme vertices) 
from all iterations of steps 6 through 9. 
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11. Repeat steps 2 through 10 for each alternative 
vertex center. 

12. Choose the local center from step 10 having the 
minimum radius as absolute center of the graph. 

F. EXAMPLES 

The following two examples will illustrate the applica- 
tion of the vertex-center-approach algorithm. Example 3 
illustrates a case in which the central path is found to 
contain a cross-edge, while Example 4 illustrates how newly- 
defined extreme vertices can prevent this. 

Figure 6, Example 3, shows a graph with the edges in Tc 
in solid lines and the cross-edge as a broken line. The ex- 
treme vertices and vertex center are appropriately marked. 
The distances between all pairs of extreme vertices, in- 
cluding non-recursive paths over cross-edge (4,5), are 



d(p,c,q) 


= d(3,4,l,5,6) 


= 35 , 


d(p,c,r) 


= d (3 , 4 , 1 , 2) 


= 24 , 


d (q , c , r) 


= d(6 , 5 , 1 , 2) 


= 23 , 


d (p , 4 , 5 , q) 


= d(3 ,4 , 5 ,6) 


= 18 , 


d(p,4 ,5 ,c,r) 


= d(3,4,5,l,2) 


= 33 , 


d(q,5 ,4 ,c,r) 


= d (6 , 5 , 4 , 1 , 2) 


= 34 . 




Figure 6: Example 3 
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It is clear from this list of path lengths that all alter- 



nate paths are shorter than (p,c,q) . Since removing either 
edge (1,4) or (1,5) breaks (p,c,q), it is ignored and only 
paths (p,4,5,c,r) and (q,5,4,c,r) need be considered, as 
outlined in Section III.D, since these are the two next 
longest paths. The longer of these two is the path (q,5,4,c,r) 
with a path length of 34, and it can be broken by removing 
edge (1,4) in the circuit (1,5, 4,1). Thus, removal of edge 
(1,4) eliminates both paths (p,c,q) and (q,5,4,c,r), and 
leaves path (p,4,5,c,r) as the longest path over the modified 
tree, with a length of 33. The midpoint of this path is the 
point Xj with a radius r(xi) = 33/2 = 16^. The reader may 
readily verify that Xi is the absolute center of the graph 
in Figure 6 by applying the Rosenthal-Smith algorithm. 

The use of relations (11) and (12) were omitted above 
for illustrative purposes. They would normally be applied 
in step 5 of the algorithm in the following manner: 



1(4,5) = 10 < d(4 , 1 , 5) = 27 ; (11) 

(12a) 

d(4 , 1) + d ( 5 ,4) - d (1 , 2) = 14 + 4 - 6 = 12) 



1(4,5) < max / > = 12 . 

ld(5,l) + d(3 , 4) - d(l , 2) = 13 + 4 - 6 = llj 

Since both relations hold for strict inequality, one knows 
that adding cross-edge (4,5) to Tc will result in a reduction 
of the radius of G providing no newly-defined extreme vertex 
generates a path longer than (p,c,q). 

Had the length of (4,5) been 12 so that equality held in 
(12a), path (p,4,5,c,r) would have been 35 units long, the 
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same as (p,c,q). In this case Xi would be an alternate ab- 
solute center, with the point x 0 , the midpoint of (p,c,q), 
being the other; and both would have a radius of 35/2 = \ 1 \. 

If 1(4,5) is increased to greater than 12, both paths 
(p,4,5,c,r) and (q,5,4,c,r) become longer than (p,c,q); 
and no reduction in the radius of G is possible. The point 
x 0 is the only absolute center in this case. 




Figure 7, Example 4, is a modification of Figure 6 with 
vertices 7 through 12 added. The distance d(l,5) is in- 
creased to 14 units from 13 and 1(4,5) is increased to 11 
units. Relations (11) and (12a) are still the applicable 
forms of the distance criteria, but now they are 
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1(4,5) 


= 11 < d(4 , 1 , 5) = 28 


(11) 


1(4,5) 


[ 14 + 4 - 6 = 12' 

< max < >= 12 . 

1 14 + 4 - 6 = 12 


(12a) 



This indicates that adding cross-edge (4,5) to Tc and break- 
ing the circuit along either path (4,1) or (5,1) will reduce 
the radius of G. Since d(q,5,4,c,r) = d(p,4,5,c,r) = 35 and 
d(p,c,q) = 36, it would seem that either (4,1) or (5,1) 
could be broken and obtain equal reductions. 

Arbitrarily choose to break the circuit by removing an 
edge in the path (4,7,8,9,c). When each of these four edges 
is removed in turn, the following distances from newly de- 
fined extreme vertices are obtained: 

• by removing (4,7), d(7,c,5,4,p) = 40 < d(p,c,q); 

" ” (7,8) , d(8,c,5,4,7) = 36 = " ; 

" " (8,9), d (9 , c , 5 , 4 , 8) = 36 = " ; 

" ” (9,c) , d(r , c , 5 , 4 , 9) = 43 < » . 

From these results it is seen that if either edge (4,7) or 
(9,c) is removed, the radius of G is increased . If (7,8) 
is removed, then d(7,4,5,xi) = 17^, where Xj is the midpoint 
of the path (p,4,5,c,r) with d(p,4,5,Xi) = d(r,c,12,xi) = 17 . 
Since the midpoint x 0 of the path (p,c,q) is colocated with 
the vertex center and has radius r(x 0 ) = r = 18, xi may be 
moved \ unit in either direction before its radius becomes 
greater than r(x 0 ). However, moving Xi ^ unit toward vertex 
5 reduces d(7,4,5,Xi) to only 19 units which is still greater 
than r(x 0 ). Similar results occur when edge (8,9) is removed 
and when one attempts to break the circuit by removing an edge 
in the path (c , 10 , 11 , 12 , 5) . 
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Thus no reduction in the radius of the graph in Fig- 
ure 6 is possible due to the generation of newly-defined 
extreme vertices, and it is concluded that the point x 0 
is the absolute center of the graph with r a = r c = 18. 

In both of the above examples the vertex center of the 
example graph has been unique. When there are alternate 
vertex centers, the algorithm is repeated for each vertex 
center as given in step 11 of the algorithm, and the opti- 
mum result(s) is(are) chosen in step 12. 

The reader should bear in mind that this algorithm will 
find the absolute center of a graph only when the vertex 
center lies on the central path of G. 

G. EXTENSION TO THE ABSOLUTE M-CENTER 

The algorithm may be extended to deal with the absolute 
m-center, m > 1 , in the same manner used by Gillespie. 

1. Find the vertex m-center, V * = {c.}, i = 1, ... ,m. 

m l 

2. Form m minimum distance trees, Tc^ rooted on c^, 

i = 1, ... ,m, by connecting each vertex in V to the nearest 

c^. If any vertex is equidistant between two or more vertex 
centers, place it in both (all) associated trees. 

3. Apply the preceding algorithm for the absolute cen- 
ter to each of the trees Tc^, i = 1, ... ,m. Do not con- 

sider any cross-edges between vertices not in the same tree. 

4. If any vertex was found to be equidistant between 
two or more of the m vertex centers in step 2, choose that 
affected tree having the smallest individual absolute radius 
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to retain the tying vertex, remove it from the other trees 
and repeat step 3 for those affected trees. 

5. Choose the maximum of the set of m individual ab- 
solute radii as the absolute m-radius r^ for this partition 
of G. 

6. Repeat steps 1 through 5 for all alternate vertex 
m-centers . 

7. Choose that vertex m-center (partition of G) yield- 
ing the minimum absolute m-radius as the final solution. 

Gillespie used the vertex 2-center as the basis for 
partitioning a graph in the first step of his algorithm, 
but avoided any consideration of its validity or uniqueness 
other than to examine alternate vertex 2-centers. The ex- 
tension of this method to the use of the vertex m-centers 
is used here with the implicit understanding that, while the 
result it yields may not be optimal, it at least gives a 
relatively easily determined upper bound for the problem 
when m is small. While no research has been done on the 
subject, it is believed that the results become less reli- 
able as m increases. No examples will be given of the ap- 
Pl ication of the above algorithm. 

Rosenthal and Smith [7] make no mention of alternate 
optima; and whereas their method of moving vertices into 
different subgraphs until no further reduction in radii is 
achieved apparently does find an optimum solution, it is 
not clear what is done when alternate optimum centers exist. 
Their method appears to minimize all individual radii to the 
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greatest extent possible rather than concentrating on the 
maximum radius, which would satisfy the definition of the 
absolute m-radius given herein. 
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IV. THE ARC-CENTER APPROACH TO THE ABSOLUTE M-CENTER 



Recall that by Theorem 2 the absolute center of graph 
G is located at the midpoint of the central path, where the 



The Rosenthal-Smith algorithm based on this definition proves 
to be an iterative process which seldom finds a non- recurs ive 
path in the first iteration. It was felt that one of the 
primary shortcomings of this method of attack is the lack of 
any definite ' relationship between the location of the absolute 
center and the middle vertex of the trio defining the central 
path. That is, the absolute center need not occur on an edge 
incident to the middle vertex. There may in fact be several 
choices for the middle vertex, all defining the same path. 

Since the absolute center of G is constrained to lie at 
some point on G, and therefore lies either at an interior 
point or end point (vertex) of some edge in G, it was decided 
that a procedure based on finding a central arc or arc-center 
(an edge containing an absolute center of G) might be more 
efficient. This led to the arc-center algorithm for the ab- 
solute center of a graph, developed in the following para- 
graphs. Since the absolute center need not be unique, the 
arc-center need not be unique. 



central path has been defined as the path (v ,v ,v ) satis- 

A V Z# 

fying 






mm 



j 
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A. FINDING THE ARC-CENTER 

Define the arc-center to be an edge in G containing an 
absolute center of G. It should be noted that, in the 
event an absolute center lies at a vertex of G and only in 
that event, there are two arc- centers (or equivalently, 
none) associated with that particular absolute center. 

Let the edge (i,j) be the arc-center of G and parti- 
tion G into two subgraphs G^(V^,A^) and Gj(Vj,Aj) such that 

all elements of V. are closer to vertex i than to vertex j 

1 J 

and all elements of V. are closer to vertex j than to ver- 
tex i as determined from the distance matrix D (See Fig- 
ure 8 ). There may be vertices which are equidistant between 
i and j and thus could be placed in either G^ or G^ . This 
event will be referred to as a type A tie, and a method of 
dealing with it will be developed later; but for the time 
being, assume that no type A ties exist. 

Let the vertices in be ordered according to distance 
from vertex k in a manner similar to that developed in Sec- 
tion III, and denote them p k> P kl >P k2 >. ••• > %> %l> “• » 

r^, ... , etc. For brevity let their corresponding dis- 

tances to vertex k be denoted d(p k ,k) = dp k , d(p k -^,k) = 



dp 



kl 



Then, associated with edge (i,j). 



dp. = d(p.,i)= maxi min 
1 1 i \ v£V 



d (v , i) , d(v, j) 



]} 



= max d(v,i) 



v cV i 



dp ,• = d(p. ,j)= max 
J J j 



( 14 ) 



mm 

veV 



d(v,i) , d(v, j) 



= max d(v,j) . 
v eVj 
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Ties between extreme vertices within and G^ may be re- 
solved arbitrarily. 

Using this notation, a nearly equivalent definition 
of the central path is that path P* in G(V,A) which yields 

1 ( P* ) = min [ dp. + dp. + l(i,j)]. (14a) 

( i > j ) e A 1 3 

If an absolute center x* occurs on any edge (i,j) then 
dp^ + d(i,x*) = dpj + d(j ,x*) ; or equivalently, d(j,x*) 

- d(i,x*) = d(p^,i) - d(pj,j). If the absolute values of 
the latter equation are taken, one obtains 

l(i>j) > I <ip ± " dp. | = |d(p i ,i) - d(pj , j ) | , (15) 

since 

l(i J) > | d(j ,x*) - d(i ,x*) | . 

On first consideration it may appear that if all M pos- 
sible partitions of G into two subgraphs (as described 
above) for each edge (i,j) in A were examined and 1(P) = 

[dp^ + dpj + l(i,j)] was computed for each partition, an 
edge which simultaneously satisfies (14a) and (15) would 
necessarily be found. The absolute center x* would then be 
the midpoint of the path P* = (jd, i,j,Pj), and it would be 
located on edge (i,j) a distance 

d(i,x*) = % [l(i,j) - dp. + dp.] (16) 

from vertex i. However, due to the existence of circuits in 
G, it frequently happens that one or more vertices are placed 
in G^ but would be closer to an absolute center if they had 
been placed in G^ . When this occurs, incorrect path lengths 
are defined. 
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For instance, in Example 6, Figure 10, the absolute 
center is colocated with vertex 1, and the shortest path 
from vertex 5 to the absolute center is along edge (1,5) 
with a length of five units. But vertex 5 is only two units 
from vertex 4 along edge (4,5), so when set (1,4) is tabu- 
lated, vertex 5 is placed in V 4 instead of Vi which makes 
d(5,l) = 7 in this partition. Furthermore, since vertex 5 
is the primary vertex p 4 in V 4 , an incorrect central path 
candidate is defined: Path (5, 4, 1,3). Had vertex 5 been 

placed in Vj, the central path candidate would have been 
the path (5,1,4) which is the actual central path of the 
graph. Note that d(5,4,l,3) = 9 < d(5,l,4) = 10. 

It is generally true that, if any vertex is incorrectly 
partitioned and would have been a primary (most extreme) ver- 
tex under correct partitioning, an incorrect central path 
candidate is defined which is shorter than the actual central 
path of the graph. A method of dealing with these cases of 
incorrect partitioning is developed below by using the ap- 
proach developed by Rosenthal and Smith in defining the cen- 
tral path. 
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Consider the representative partition of G into {G^,G^.} 

shown in Figure 8 and note that there are three possible 

candidates for the central path shown in the figure. These 

are the paths (p-^ij.p.), (p^i^), and (PjJ.qj)- If 

the path (p^,i,j,p.) is the actual central path of G, one 

of two possible conditions will exist; either the absolute 

center is on (i,j) (in which case (15) will be satisfied by 

(i,j)) or the absolute center is on another edge such as 

(b,i) or (j ,d) (in which case (i,j) cannot satisfy (15)). 

In either case, however, the path (p^,i,j,Pj) must satisfy 

(14a) . In the latter case it is impossible to predict what 

conditions will arise when G is partitioned into (G, ,G.} 

b i 

and {G.,G.} without knowing what the connective pattern and 
J a 

distances are in G, but it is evident from Figure 8 that the 

set G as it occurs in {G.,G.} differs from the set G. that 

i J i 

would occur in {G^,G^}. It will become apparent in later 
examples that this is generally true. 

Now assume that d(p^,i,q^) > d(p^,i,j,p^) in Figure 8. 
This is equivalent to 



d(q i ,i) > d (p ^ , j ) + 1 (i , j ) . (17) 

When this occurs, edge (i,j) cannot satisfy (15) because (17) 
must also hold if d(p^,i) is substituted for d(q^,i) in (17); 
but this requires the midpoint of the path (p^,i,j,Pj) to lie 
on the branch (p^,i), not on the edge (i,j). When the dis- 
tances dp^, dq^, and dp^. for a given partition of G are such 
that (17) holds, this will be referred to as a ’B condition’. 
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Thus, when a B condition occurs, the path (p^,i,q^) is the 
longest non- recursive path through vertex i, and so be- 
comes a primary candidate for the central path of G as de- 
fined by Rosenthal and Smith [ 7 ]. Therefore, if (p^,i,q^) 
is the shortest path found, it is the cental path of G, and 
the absolute center of G is the midpoint of (p^,i,q^) with 
absolute radius 



r a = 1 5[d(p i ,i) + d (q i , i) ] . 

Vertex a in Figure 8 is intended to be equidistant be- 
tween vertices i and j; a type A tie. Under the condition 
of the figure this creates no problem, but had vertex a 
been either p^ or p^ the paths (a,i,j,Pj) with a in G^ and 
(P a »i,j»a) with a in G^ (assuming a can be either p^ or p^) 
would probably be of different lengths. Should this occur, 
both paths must be checked and the shorter chosen as the 
central path candidate; since by definition, if the central 
path is either of these paths, it must be the shorter of the 
two . 

It is also possible that (i,j) does satisfy (15) when 

vertex a is in G. , but does not when a is in G . ; or that a 

1 3 

enters into a B condition as p^ or q^ when in G^ ; but (15) 
is satisfied when a is in Gj , etc. The algorithm checks 
all possible combinations of type A ties and B conditions, 
and selects the shortest central path candidate as the cen- 
tral path of G. 
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B. THE ALGORITHM 

Using the results of the preceding discussion, the arc- 
center algorithm for the absolute center of G may be formu- 
lated as follows. 

1. Find the minimum of d(v,i) and d(v,j) for all M 
pairs of adjacent vertices by working from the distance 
matrix D using the same technique as for the vertex 2-cen- 
ter. This gives M different partitions |g^ ,G^ of G. When 
doing hand calculations, the results are best tabulated in 
an (n+5) x (M+l) tableau as follows (see later examples, 
also) . 

a. List all edges (pairs of adjacent vertices) of 
G in row 1, c.olumns 2 through (M+l); i.e., across the top of 
the tableau. 

b. List all vertices in G in column 1, rows 2 
through (n+1) . 

c. Place the following labels in column 1 rows 

(n+2) through (n+5) respectively: ’ dp^ ,dp^ ’ , *1 (i , j ) ' , 

* 1 (P) f , 'B'. 

d. In row (n+3) (labeled ’l(i,j)') write the length 
of each edge in the column headed by the pair of adjacent 
vertices defining that edge. 

e. Divide the column under each pair of vertices 
into two sub-columns in rows 2 through (n+1) . Apply the 
vertex 2-center algorithm to G for each pair of vertices in 
row 1. Enter the minimum value of [d(v,i), d(v,j)] for each 
vertex v in the row labeled with that vertex and the sub- 
column headed by the vertex i or j for which the minimum occurs 
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(the cell of the tableau corresponding to v and i or j, 
whichever minimizes [d(v,i), d(v,j)]). If d(v,i) = d(v,j), 
enter this value in both subcolumn i and subcolumn j and 
mark this as defining, a type A tie by placing 'A' in a con- 
venient location, such as between the two entries. 

f. When completed with step e, examine the mini- 
mum distance entries in each set (column) and, if the maxi- 
mum entry in one subcolumn is less than or equal to the second 
or higher ordered (lesser valued) maximum entry in the other 
subcolumn of the set, mark the set as having a possible B 
condition by placing * b T in a convenient location, such as 
beside the smallest entry which is greater than or equal to 
the maximum entry in the other subcolumn. 

2. For each partition (set or column) find the max 

veV^ 

d(v,i) = dp. and the max d(v,j) = dp.. 

1 veV . J 

1 

a. If a set has been marked as having a type A tie 

and the tied value is either dp^ or dp^ , delete the tying 

value from and find the resulting values of dp^ and dp^ , 
then move the tying value from to Gj and recompute dp^ 
and dp.. This defines two different sets {G.,G.} wherein 
the vertex having the tied value is in G^ for one set, and 
in for the other. 

b. Should there be more than a single pair of type 
A ties, say k pairs, in a given set, each of the 2k possible 
significant combinations is tried as follows: 

1) The maximum tying value is deleted from G^ ; 

all others remain in G^. Check that the resulting route is 
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a non-recursive path and if it is, complete the calcula- 



tions. If it is a recursive path, the tied values now in 
are moved to sequentially, and the route is checked 
for a non-recursive path each time; only the combination 
which determines the longest path is retained, all others 
being dropped from consideration since this is part of the 
maximizing step in finding the central path. 

2) When step 1) is completed, move the maxi- 
mum tied value to G^ and repeat step 1) , moving the others 
to G^ sequentially. When this process is completed, select 
the non-recursive path(s) with the maximum length dp^ + dp^ 
+ l(i,j), including ties, and use it (them) as the path(s) 
representing the set(s). Ties yield alternate paths and, 
in the event their path lengths are the minimum of all cen- 
tral path candidates, they will yield alternate absolute 
centers . 



3. For each partition {G^ ,Gj } compute 1(P) = dp^ 

+ dpj + l(i,j) to find the total path length. 

4. Choose the minimum value in the set of total path 
lengths. If this is greater than the length of any legiti- 
mate central path candidate previously resolved, go to step 
8 and terminate. If not, test the set with the distance 
criterion (15). 

5. a. Go to step 6 with any set(s) (including ties) 
which satisfy (15) . 

b. If the minimum length path found in step 4 
fails to satisfy (15), check for a B condition. In filling 
in the tableau possible B conditions are marked; i.e. , if 
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the second maximum entry is subcolumn i is greater than the 
maximum entry in subcolumn j , the set is marked as having a 
possible B condition. Let dki be the maximum entry in sub- 
column k, dk 2 the second maximum entry, etc., k = i, j. 

Then if di 2 > dj j + l(i,j) it is possible that a B con- 
dition exists, but it is necessary to check that a non-re- 
cursive path exists from the vertex yielding dij to vertex 
i and then to the vertex yielding di 2 . Note that dki is al- 
ways dp^, k = i, j (in step 2.b. only the maximum value was 
considered in each subcolumn, so it was not necessary to 
consider this problem at that time). When the value of dk^ , 
h = 2, 3, ... , has been found which yields the longest non- 

recursive path (p^.,k,q^.) (where q^. is the vertex yielding 
this value of dk^) , check that this dk^ still satisfies (17) 
i.e., for k = i, check that di^ = dq^ = d(q^,i) > d(pj,j) 

+ 1 (i > j ) • 

1) If this holds, a B condition exists; the 
local radius is Jg(dp^, + dq-j.) and the local center is located 
at the midpoint of the path (p^,k,q^.) where k refers to the 
subset or Gj containing the B condition. Enter the sum 
(dp^ + dq^) i n the ' B ’ row of the tableau. Return to step 4 
and choose the next minimum total path length. If this is 
longer than the B condition path (dp^. + dq^.) just found, go 
to step 8 with the results from the B condition set just re- 
solved and terminate. 

2) If dq^ < dpj + 1 ( i , j ) , discard the set; re- 
turn to step 4 and choose the set having the next minimum 
1(P). 
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6. Divide the total minimum path length by two to find 
the local radius. 

7. If no condition B was involved, locate the local 
center x* at a distance d(i,x*) units from vertex i along 
(i,j) where d(i,x*) is as defined in equation (16). If a 

B condition was involved, the local center is found as given 
in step 5 .b . 1) . 

8. The algorithm terminates when the set with minimum 
total path length satisfies (15) , or when it is verified that 
the condition B path of minimum length is found and is shorter 
than or equal to the shortest path satisfying (15) . 

C. EXAMPLES 

Hakimi's example (Figure 1) will be used first to il- 
lustrate the application of the algorithm. The graph in 
Figure 1 is reproduced in Figure 9 for convenience, along 
with the associated distance matrix and the computational 
tableau resulting from applying the arc-center algorithm to 
locate the absolute center of this graph. 

Upon completion of steps 1, 2, and 3 of the algorithm, 
the tableau appears as in Figure 9 with the exception that 
the 'B' row is not filled in. The only type A tie involved 
occurs .in set (2,3) for vertex 4. Step 2. a. resolves the 
tie by first placing vertex 4 in G 2 (deleting the entry *12' 
from subcolumn 3 of set (2,3) which results in the path 
(4, 2, 3, 5) with dp 2 = 12, dp 3 = 10 and 1(P) = 36; vertex 4 is 
then moved to G 3 (the entry '12' is moved to subcolumn 2) 
which results in the path (1,2, 3, 4) with dp 2 = 10, dp 3 = 12 
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Figure 9: Example 5. 
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and 1(P) = 36. The entry '36' is made twice in row ' 1(P)', 
once for each path found in step 2. a. 

The minimum path length (minimum entry in row ' 1(P)') 
found in step 4 is 34, which occurs for sets (1,2) and (3,4). 
Step 4 selects one of these sets, say set (1,2) and finds 
that it fails the distance criterion (15); i.e., 1(1,2) = 

10 < | dp x - dp 2 | = |0 - 2 4 1 = 24. Step 5.b. notes that set 
(1,2) is marked as having a possible B condition and finds 
that dk 2 = d2 2 occurs for vertex 3 with an entry of 14, but 
the path (3,2,5) is recursive (paths (3,2) and (2,5) have 
the edge (2,3) in common), thus vertex 3 cannot be the q2 
vertex even though it does satisfy (17) . The third maximum 
entry in subcolumn 2 occurs for vertex 4 with dk 3 = d2 3 = 12. 
The path (5,2,4) is non- recursive ; therefore vertex 4 is the 
q 2 vertex. Furthermore, inequality (17) still holds, so a 
B condition does exist. 

Step 5.b.l) then finds the local radius to be ^(dp 2 
+ dq 2 ) = Jg(24 + 12) = 18. The center of the path (5,2,4) is 
found to occur at point y on the graph in Figure 9, located 
6 units from vertex 3 on edge (3,2) . The value of dp 2 + dq 2 
= 36 is entered in the * B * row of column (1,2) and the algo- 
rithm returns to step 4 to select the next minimum path length. 

The next minimum path length is also 34 and occurs for 
set (3,5) as previously noted. Conditions very similar to 
those just discussed for set (1,2) are found to exist in set 
(3,5), with vertex 4 being the q vertex again. In this case, 
however, the B condition path is found to be the path (4,3,1) 
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with midpoint at the point x on edge (2,3) of the graph; 
and the radius is again found to be 18. 

Returning to step 4, the next minimum value of 1(P) = 

36 is found to occur for both partitions under set (2,3). 

The two paths defined are (4, 2, 3, 5) and (1,2, 3, 4), as pre- 
viously noted; but these are the same paths defined by the 
B conditions in sets (1,2) and (3,5) respectively, so set (2,3) 
may be dropped from consideration. If they were not dropped 
they would be found to satisfy (15) , and thus they give the 
two alternate central paths and absolute centers directly 
without any consideration of B conditions. 

The next sets eligible for consideration are (2,4) and 

(3.4) with 1(P) = 42 for each, but this is strictly greater 
than the lengths of the B condition paths previously re- 
solved for sets (1,2) and (3,5), so the algorithm terminates, 
having found the alternate central paths (4, 2, 3, 5) and 

(1.2. 3. 4) with midpoints y and x respectively. The absolute 
radius is then r(x) = r(y) = r = 18. 

Example 6, Figure 10, illustrates a case in which the 
central path cannot be defined by (14a) due to the existence 
of a circuit, as previously discussed. This case is charac- 
terized by a B condition in the minimum-path set as well as 
by the set's failure to satisfy (15). 

The first iteration chooses either set (1,2) or (2,3) 
as the minimum-path set. If set (2,3) is chosen first, it 
is found to not contain a valid B condition as the paths 

(2.1.4) and (2,1,5) have the edge (1,2) in common. The set 
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Figure 10: Example 6. 
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is therefore dropped from consideration and set (1,2) is 
chosen on the second iteration. 

Set (1,2) defines the same initial paths as does set 
(2,3), but now the possible B condition is found to yield 
the valid B condition path (4,1,5) with a length of 10 
units. The center of this path lies at vertex 1 and has a 
radius of five units. 

The third iteration finds in step 4 that all remaining 
sets have 1(P) > 10, and so the algorithm terminates in 
step 8 with an absolute radius of 5 and the absolute center 
at vertex 1. 

Note in Figure 10 that the set (4,5) contains three 
type A ties ,- associated with Vertices 1, 2, and 3. When 
step 2.b. is called on to resolve these, the paths (3, 2, 1,4, 5) 
and (3, 2, 1,5, 4) are found. These two paths occur when all 
type A tied values are in the same subcolumn, and are the 
same paths defined by sets (1,4) and (1,5) respectively. 

Example 7, Figure 11, illustrates a less obvious case 
generating both a type A tie and a B condition. Set (1,2) 
has a type A tie corresponding to vertex 3 and a B condition 
(which is never fully resolved) for the partition with ver- 
tex 3 in Vi (subcolumn 1) of set (1,2). The first iteration 
through step 4 chooses set (3,4) which is found to have a B 
condition involving vertices 1 and 2 which yields the central 
path candidate (1,3,2) with a path length of 10. This gives 
a local center at vertex 3 with a corresponding radius of 5. 

The second iteration through step 4 finds the partition 
in set (1,2) with 1(P) = 10, but it is marked as having a 
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Figure 11: Example 7. 
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possible B condition. When this is resolved the B condi- 
tion path (4,1,3) is found with a path length of 11, but 
this is strictly greater than the length of the central 
path candidate found in the previous iteration, so the set 
is dropped from further consideration. 

The next iteration through step 4 shows that all re- 
maining sets have 1(P) values greater than the path found 
in the first iteration (10); therefore, the algorithm ter- 
minates in step 8. 

The Rosenthal-Smith algorithm encounters trouble in 
dealing with peripheral vertices, where a peripheral vertex 
is defined as any vertex which is connected to the remainder 
of the graph via only one other vertex of the graph. In 
the calculation of the central path (v ,v ,v ) , if v is a 
peripheral vertex, all central path candidates associated 
with it have backpaths. The following three examples will 
illustrate how the arc-center approach eliminates this pro- 
blem. 

In Example 8, Figure 12, vertex 4 is a peripheral ver- 
tex and is connected to the rest of the graph by edge (2,4) 
via vertex 2. The arc-center computation tableau entries 
are all zero in subcolumn 2 of set (2,4) as they must be for 
any peripheral vertex. 

The first iteration of the algorithm through step 4 
chooses set (2,4) with a path length of 5 and resolves the 
possible B condition to find the B condition path (1,2,3) 
with a length of 6 units. This gives the local center x* 
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Figure 12: Example 8. 
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on edge (1,2) of the graph in Figure 12, with a local radius 
of 3 units. 

The three remaining sets all have 1(P) = 6, so all must 
be resolved. Set (1,2) is found to not contain a valid B 
condition, and so defines the path (3,2,1). The set satis- 
fies (15) , so it is retained and found to again yield the 
local center x* . 

The algorithm finds that set (1,3) does not contain a 
valid B condition but does satisfy (15) also. In this case 
the set is found to define the path (4, 2, 3,1) with midpoint 
y* located at vertex 3 and a radius of 3 units. 

Set (2,3) is found to satisfy (15) on the last itera- 
tion. The path defined is (4, 2, 3,1), the same as that de- 
fined by set (1,3). Therefore, the absolute radius of the 
graph is 3, with alternate absolute centers, points x* and 
y*. 

In this example, 3 arc-centers were found. The first, 
edge (1,2), was found through a B condition path first, and 
then through its own set, which satisfied (15). Since the 
alternate absolute center y* occurs at a vertex, two arc- 
centers, edges (1,3) and (3,2), were found for it. 

Example 9, Figure 13, is a variation of Example 8 with 
1(2,4) increased to seven units. In this case sets (1,2) 
and (2,4) tie for minimum path length with 1(P) = 11, but 
set (1,2) does not satisfy (15) while set (2,4) does. Both 
sets define the path (1,2,4) which has its midpoint at point 
z*. The remaining sets both have 1(P) > 11, and so the 
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Figure 13: Example 9 



algorithm terminates with z* as the absolute center and 
r & = Sh. There were no B condition sets. 

Example 10, Figure 14, illustrates a case in which a 
vertex (vertex 4) is technically not peripheral because it 
is connected to two other vertices, 1 and 3, in the graph. 
However, 1(1,4) is greater than 1(1,3) + 1(3,4) (8 > 2 + 2 
= 4) , therefore all minimum distance paths to vertex 4 must 
pass through vertex 3, and vertex 4 behaves as a peripheral 
vertex. 

There are three sets which tie with a minimum path 
length of six. These are: (1,3) and (3,4) which fail to 

satisfy (15) but do not contain B conditions; and (2,3), 
which does satisfy (15). Set (2,3) defines two alternate 
paths, (2,3,1) and (2,3,4), both having the same path length; 
but since (2,3) satisfies (15), the absolute center is on 
edge (2,3) and it is immaterial which path is considered. 

The absolute center is therefore a point w and the absolute 
radius is 3. 

Examples 9 and 10 illustrated that no further special 
consideration need by given to peripheral vertices. This 
is a definite improvement over the Rosenthal - Smi th algorithm, 
which, as previously noted, requires special inputs for 
peripheral vertices when programmed into a computer. 

D. EXTENSION TO THE ABSOLUTE M-CENTER 

Define the arc-m-center of the graph G(V,A) to be a set 
of m edges A m * such that for every other set of m edges A^eG, 
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Figure 14: Example 10. 
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( 18 ) 



max d(v,A ) > max d(v,A *) = r ’ . 

,, m - m m 

veV veV 

Where d(v,A m ) is the minimum distance from vertex v to any 

endpoint (defined vertex) of the set A m » Since every edge 

has two endpoints, A^ defines a set of m pairs of adjacent 

vertices V^. Inequality (18) can therefore be written 

in a form consistent with (5) : 

max d(v,V ? ) > max d(v,V *) = r ’ . (19) 

veV veV z 

1 (A m *) may now be defined by generalizing (14a) in the fol- 
lowing form (using the same notation and definitions) : 

1(A*) = min f max [dp. + dp. + l(i,j)]y. (20) 

A|n £ A \ J J 

Equation (20) outlines an iterative algorithm for par- 
tioning G into a set of m subgraphs {G^}, i = 1, ... ,m, 

and simultaneously finding the central path and hence the 
absolute radius r(G^) , of each G^. As the algorithm iterates 
through all possible partitionings of G (A A) , that parti- 
tion yielding the minimum absolute radius is found. Since 
the set of points X m * yielding the minimum absolute radius 
defines an optimum partition of G {G^*}, and since the algo- 
rithm will test all possible partitions of G, the algorithm 
must find {G^*}. 

The computational procedure is essentially the same as 

previously outlined for the single arc-center algorithm, and 

M 

iterates through ( m ) modified arc-center calculations. The 
optimum result from the set of all results is then chosen. 

The changes required in the arc-center algorithm are 
listed below. 
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1. In step 1, the algorithm now starts by choosing 

one of the f^) combinations of pairs of adjacent vertices 
m 

not previously examined; these are- the sets V ^ • 

a. Since there are k < 2m different vertices in- 
volved, the computation tableau may be reduced by k rows. 

This is possible because all entries in rows corresponding 

to elements of the set under consideration have zero entries. 

b. The minimum distance from each vertex (row in 

D) to any element of is entered in all subcolumns cor- 

responding to that element. This produces ties any time 
two or more edges sharing a common vertex are in A . Call 
these a type C tie. 

2. No Change. 

a. No Change; but there could be a type A multi- 
ple tie; i.e., one vertex could be equidistant between 
several elements of • It is still allowed in only one 
set at a time, and checked in the same manner as in the 
arc-center algorithm. 

b. No Change; same comment. 

3. Compute 1(P) for each subset with and without type 
C ties. For each group of sets with a common type C, choose 
one with, rest without to min(max 1(P)). 

4. Examine all tying sets generated by procedure A in 
step 2 and choose those with minimum path length for further 
consideration. Do not discard the others at this time. 

5. a. No Change. 

b. If no condition B exists, retain the subset as 
originally defined, and process it as though it had satisfied 
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(15) except the center of the path is simply the midpoint 



of (P^i > j >Pj) • 



X 




J 






1) . No Change . 






2) . Delete . 


6. 


No 


Change . 


7. 


No 


Change . 


8. 


If 


all m subsets have not been processed, return t 


step 4. 






9. 


If 


any type C ties remain, discard the entire set 



and return to step 1. 

10. Find the maximum radius in the set of m radii. 

a. If this radius is less than the maximum radius 
retained from any previous iteration, discard all previous 
results and store all results of this iteration. 

b. If this radius equals that retained from any 
previous iterations, retain the previous results and store 
the present results also. 

c. If this radius is greater than that retained 
from any previous iterations, discard the results from this 
iteration . 

11. If all (M) combinations have not been examined, re 

m 

turn to step 1. 

The analyst now may choose between any sets of m-cen- 
ters which may have tied for min(max radius) , using some 
further criterion. 

E. EXAMPLE 

Example 11, Figure 15a, is an abbreviated absolute 2- 

center problem which will illustrate the application of the 
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arc-m-center algorithm. Only six of the possible (^) = 66 
sets ^ 2*2 are s ^ own i- n Figure 15b. Note that in each set 
(iteration) , the pair identity of the vertices is retained; 
they cannot be mixed randomly. Since three sets are shown 
side by side in each tableau, it was not possible to delete 
as many rows from the tableaux as would normally be done if 
a separate tableau were made for each set. 

Observe that in set 1, the set {(1,2), (4,6)}, vertex 
7 is equidistant between vertices 1 and 6. While this is 
not a type C tie as previously defined since two different 
vertices are involved, it is handled in the same manner. 

In this instance subset (4,6) is not affected when vertex 7 
is removed; and since its path length is greater than that 
of subset (1,2) with vertex 7 included, it makes no dif- 
ference which subset includes vertex 7. 

Set 5 is entirely different. Here vertex 4 is equi- 
distant between vertices 5 and 6, producing both type A 
and C ties. When vertex 4 is in subset (6,7), the path 
(4,6, 7,1) is defined with a length of 11; and (5,6) de- 
fines the type B path (2,5,3), with a length of 13. When 
vertex 4 is moved to subset (5,6), a type A tie is pro- 
duced which is resolved into path (4, 6, 5, 2) with length 
16, which satisfies (15); or the above mentioned type B 
path (2,5,3), with vertex 4 being "peripheral" to the path. 

Set 4 is the set which satisfies min(max radius) of 
those sets shown in Figure 15b and may be shown to be the 
optimum set. Subset (1,2) defines G! = {1,2,7} with abso- 
lute center at point x; and (3,4) defines G 2 ={ 3 , 4 , 5 ,6} with 
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Figure 15a: Example 11. 
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Figure 15b: Example 11 
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absolute center at point y*, colocated with vertex 4 
have an absolute radius of 5. 



Both 
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VI. CONCLUSION 



A . SUMMARY 

A brief review of known previous algorithms for the 
solution of the absolute m-center problem has been given, 
and mention made of the problems encountered in each. In 
addition, two new algorithms were presented. These algo- 
rithms are applicable to the solution of optimal location 
problems which can be structured in a graph theoretic man- 
ner, and so posed as absolute center problems. 

The first of the algorithms, the vertex-m- center ap- 
proach, was presented only as an aid in understanding the 
problem and for use in hand calculation of solutions to sim- 
ple examples, particularly with m = 1. It is probably less 
efficient than the Rosenthal-Smith algorithm, but it does 
provide a basis for understanding and judging other algo- 
rithms . 

The second algorithm presented was the arc-m-center 
algorithm. It was originally formulated as a means of by- 
passing the path check requirements of the Rosenthal-Smith 
algorithm. As the reader has seen, the effort was unsuc- 
cessful to the point that, for m = 1, the algorithm is quite 
Possibly significantly less efficient than the Rosenthal- 
Smith algorithm. However, for m > 1, it is felt that the 
arc-m-center approach may prove more efficient and more 
easily implemented than any previous algorithm. This algo- 
rithm can be made more efficient, both in terms of data 
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management and of required operational steps, than outlined 
herein (see following section) ; and a final choice between 
it and the Rosenthal -Smith algorithm, or a meld of the two, 
would necessarily be based on operational results. 

B. AREAS FOR FURTHER CONSIDERATION 

As previously noted, the vertex-m- center approach was 
extended to cases with m > 1 by an unproven method. Further 
research is needed to determine under what conditions the 
use of the vertex-m-center to partition a graph may yield 
an optimum partition, and a more efficient means of select- 
ing the optimum partition when the vertex m-center is not 
unique . 

This method does provide an upper bound on the solu- 
tion which might be useful in conjunction with some other 
algorithms to eliminate some of the possible solutions as 
being infeasible, and thus shorten the number of computa- 
tions and time required. It might also provide a more ef- 
ficient starting point for the Rosenthal-Smith algorithm 
than the m-node divisional path presently used. There is 
certainly a need for more investigation into these areas. 

Another approach to solving the absolute m-center pro- 
blem which might be worthy of investigation is to consider 
each of the (n) partitions of G about the central vertices 
of the m central paths, using the vertex m-center radius as 
an upper bound to reduce the number of iterations which are 
carried to completion. 
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In all previous work on the absolute m-center problem 
except that of Rosenthal and Smith (see comment, Section 
III.G), the only criterion has been to minimize the maxi- 
mum value of the set of m radii. This frequently yields a 
number of alternate solutions, all having the same maximum 
radius with different lesser radii. It is felt that atten- 
tion should be given to further criteria for choosing the 
optimum solution. Two possible considerations are to mini- 
mize the sum of the radii or to minimize the absolute dif- 
ferences between radii. 
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